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I. INTRODUCTION 

Recent astronomical measurements of Super Novae la 
light curves versus their red shifts (SNIa) Cosmic 
Microwave Background Radiation anisotropy (CMBR) 
by Wilkinson Microwave Anisotropy Project (WMAP) 
[4!|, inhomogeneous correlations of baryonic matter by 
Sloan Digital Sky Survey (SDSS) and 2dF Galaxy Red- 
shift Survey with a high precision enforce the following 
picture of cosmology: 

• the Universe is flat, 

• its evolution is consistently driven by cosmological 
constant A and cold dark matter (CDM), that con- 
stitutes the ACDM model. 

Irrespective of dynamical nature for such substances, 
that could be different, at present any model of cosmol- 
ogy has to demonstrate its tiny deviation from the ACDM 
evolution at late times, i.e. the behavior of Hubble con- 
stant should scale extremely close to 



= Hi 



(1) 



with Ho = H{to) denoting the present day Hubble con- 
stant, a — a{t) is the scale factor in the Friedmann- 
Robertson- Walker metric 



ds2 ^ dt^ - a^{t) [dr^ 



sin'^ Odip'^ 



(2) 



conveniently normalized by a (to) = 1, so that H — a/a 
with dot meaning the differentiation with respect to time 
t. The fractions J^a and JIm represent the cosmological 
term and pressureless matter including both baryons and 
cold dark matter. In ([1]) we neglect contributions by 
radiation fractions given by photons and neutrinos. 

Dynamical models closely fitting the above behavior of 
Hubble constant include a quintessence , a scalar filed 
(j) with slowly changing potential energy V{(j)) imitating 
the contribution of cosmological constant (see recent re- 
view on the reconstruction of dark energy dynamics in 
[3| ) • In present paper we find a potential of scalar field 
which exactly reproduces the scaling behavior of Hubble 



constant in flat universe with cosmological term^ (Section 
II). The potential is the square of hypersine with some 
tuned values of normalization and slope. In Section III 
we study the stability of scaling behavior versus the pa- 
rameters of potential in terms of autonomous system of 
differential equations possessing critical points. We find 
the late time attractors that further evolve to future at- 
tractors generically different from those of late time. We 
discuss a physical meaning of attractors in Section IV. In 
Conclusion we summarize our results. 



II. EXACT SOLUTION WITH SCALING 
BEHAVIOR 

The evolution is described by following equations: 



+ 3i/0-H^ = O, 



(3) 



where ps is the energy density of baryotropic matter 
with pressure pb = wbPb, which satisfies the energy- 
momentum conservation 



Pb + 3H{pB +Pb) = 0, 
yielding the scaling behavior 



PB = 



qZ(1+wb) ' 



(4) 



(5) 



We suppose the following scale dependence of Hubble 
constant 



H' = H^ 



,3(i+il.b) y ' 



(6) 



^ The other approach of reconstructing an effective potential of 
variable cosmological term by the given luminosity distance or 
inhomogeneity growth factor was considered by A.Starobinsky 
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where fls denotes the present day fraction of substance 
composed of baryotropic matter, cold dark matter and 
quintessence, that could simulate dark substances ordi- 
nary introduced in the standard consideration: the dark 
energy and dark matter. At late times (or at present) 
wb = corresponds to nonrelativisting matter with neg- 
ligibly small pressure (the dust), while wb = 1/3 stands 
for the radiation era of hot matter. 
The critical density pc is defined by 



■Pc 



(7) 



so that the baryonic and dark matter fractions are 
presently given by 



Pc 



DM 



Pdm 



Pc 



while for brevity we put 
^B — + ^DM , 



or pg ^ Pj^+ Pjj^j. 



The scalar field density and pressure 



1 



V, 



- V. 



determine the fraction 

^4, = — P0 

Pc 



1 



and state parameter function 

J 

= - 



(8) 



(9) 



(10) 



(11) 



(12) 



The substance fraction is the sum of matter and field 
fractions 



We define the vacuum energy by 

Vb = Pc^^A, 

so that in flat universe 

flA + ^S ^ 1- 



(13) 
(14) 
(15) 



Since we investigate the scaling behavior of functions 
homogeneous with respect to scale factor a, it is conve- 
nient to introduce the following variable 



N = Ina(t), 



(16) 



so that the differentiation with respect to time denoted 
by dot is reduced to the differentiation with respect to N 
denoted by prime. 



dV V 



(17) 



Then, the equation of motion for is deduced to 

^{{cp'Hry + 3{<P'Hr + r ^0. (18) 
The consideration suggests that the potential scales as^ 



V = Vo + 



fj3(t«B+l) 



Pc, 



(19) 



where fi^ is a constant. Therefore, the derivative of po- 
tential scales, too, 

V = ^3iwB + l)iV- Va). (20) 

According to (flS]) . we suggest 

3{(f>'Hy = -cV, (21) 

where the constant c can be found from ()18|) . since 

{{^'Hfy = -^V" ^c{wB + l)V', (22) 

hence, yields 
1 



c{wB + l)-c+l}V' ^0 



that is satisfied at 



One can easily get 



1-wb 



-{<t>) ^-{<t>H) 



so that 



that gives the relation 



Orf, — 



j3(«;b + 1) I 



Wb 



1 - Wb 



(23) 



(24) 



(25) 



(26) 



(27) 



In order to restore the potential yielding the scaling 
behavior, we have to resolve (|2ip making use of (|19p . 
(El, i.e. 



, 2 2pcf^0 l + WB 1 

' ) ^ 



1-WB a^^'^B+i) ' 



(28) 



^ Symmetries of evolution equations were systematically investi- 
gated in wherein the authors found the analogous behavior 
of energy density versus the scale factor and used it for deriving 
the equation of state on the basis of supernovae la data, while 
our goal is the potential itself. 
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where 



= — Pc{^K + 



1 



The integration straightforwardly yields 

^ « (0 - 0.) = arcsinh^IIlir^, 
where = SttG, and 



A = W 3(1 + wb. 



(29) 



(30) 



For brevity of formulae we put the integration constant 
(j>i, = with no lose of generality. Then, 



where 



while 



V = Voil + 



^4 



sinh 




ns 



3 1 - uil 
2 A2 ■' 



H"^ = i?o cosh^ i '2 ^ ' 



(31) 



(32) 



(33) 



So, at the present day we have cosh^ K(/>o} = 

Summarizing the result, we emphasize that there is the 
exact solution for the scalar field potential PT|) . which re- 
produce the scaling behavior of Hubble constant in the 
evolution of flat universe in presence of cosmological con- 
stant. 

The form of potential differs from the case of zero cos- 
mological constant, where the potential is the exponent 
as was studied for the scalar field with the standard ki- 
netic term in [lOj [HI [1.2, 13], while the consideration for 
the general scalar field was developed in [ij, UBI . One 
can easily notice that the present derivation is consistent 
with results concerning for the case of zero cosmological 
constant. Indeed, the integration of (|28p with the Hubble 
rate at f^A = straightforwardly gives the field propor- 
tional to the logarithm of scale factor, cj) ^ In a, that 
makes the scaling behavior of potential V cc exp{XK(f)}. 

The potential derived can be represented in the form 



V = Vo + ^Vo{cosh{XK<l)} 



1 



(34) 



at Vo = 3Vo{l — w^)/2y. So, function is composed 
by the sum of two exponential potentials with the oppo- 
site slopes and the constant positive shift of minimum. 
Such kind of potentials was investigated recently. 

In review [l6| authors presented the exact solution for 
constant parameter of state for the dark energy in the 



presence of dust-like dark matter, but not adding the 
cosmological constant, which is imitated by the dark en- 
ergy instead. In [l3| the evolution of scalar field with 
the hyper-cosine potential was studied in the case of zero 
cosmological constant: the exponential form found to be 
dominant at early times, while the square term did signif- 
icant at late times. It is clear that the late time dynamics 
essentially changed by the presence of cosmological term, 
of course. 

Various aspects of cosmological picture due to a po- 
tential given by a power of hyper-sine was investigated 
in [l8[ . The tracker properties of such the potentials were 
stressed. 

Another approach was presented in [l^ [l^, HH, [2^ . 
where authors fixed the scaling behavior of Hubble con- 
stant to find exact solutions for the scale factor a{t) in 
order to study characteristics of universe evolution. The 
consideration of recovers the scale factor behavior 
in the case of ACDM, however, the authors did not ad- 
dress the question on the scalar field potential reproduc- 
ing such the scaling. This question was investigated in 
[20], where the potential with the same form of ([M)) was 
deduced in a particular case of VIs/^a = 1/ sinh^ 1 and 
WB = 0. At this choice f^A = sinh^ l/(l-|-sinh^ 1) « 0.58, 
which is in contradiction with the recent measurements 
[j] yielding Oa — 0.766 ± 0.035. A cosmological ex- 
ploration of potential composed by sum of two expo- 
nents with opposite slopes but generically different nor- 
malization factors and a negative shift of minimum was 
investigated in [2l|, by the same method of exact 
time dependence. The authors found an oscillation of 
a{t) around some scaling dependence with oscillating 
within ]-l; -|-1]. 

In [2J] the sum of two exponents with identical nor- 
malization factors but slopes, which can be different, was 
considered at first. The late time behavior of scalar field 
energy scales both the radiation and dust, while near the 
present and future the state parameter has decaying 
vibrations around —1. Such picture differs from that of 
• Questions are the foUowings: i) What is a reason 
for the difference? ii) What can we say about a stabil- 
ity of late time and future scaling? Hi) Does presented 
exact scaling solution corresponds to fine tuned values of 
normalization and slope? These questions were not in- 
vestigated in references mentioned. We address them in 
Section III. 



III. ATTRACTORS 

Let us consider the evolution of flat universe in pres- 
ence of scalar field with potential 



V = Vo + Vo sinh^ <! ^ K, 



(35) 



where Vq, Vq and A are free parameters, which are not 
fixed by values in (|5D|) . ([?!]) . For definiteness we put 
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all parameters to be positive: Vq > 0, Vq > 0, A > 0, 
while the consideration for cases of negative values can 
be rather straightforwardly obtained from the formulae 
below. The Hubble constant is given by 



H' = — l^PB + -{^r + Vo + Vosmh'i^-Kq 



so that we introduce quantaties Uq and [/, so that 

,2 



TP — 



K 



(36) 



Then, the phase space of system is described by dimen- 
sionless variables 



^ = ^~rr^ FT ' ■^ = 77' (37) 



^/6 C/ ' ^ V3 U 
while for convenience we introduce 

1 dV V2 



U 



H 



(38) 



This choice of variables follows the observation of scaling 
in previous section: the kinetic energy and potential each 
scale like the Hubble constant squared after the subtrac- 
tion of term caused by the cosmlogical constant, while 
the derivative of potential with respect to the field scales 
like the Hubble constant squared itself. 
The definition of implies 



2,2 1 1^'^ PB 



which yields the constraint 



(39) 



(40) 



The dynamical state parameter of field is determined by 



P4> + Vq _ 



y 



Pel, - Vb + J/2 

In addition, the equations of motion produce relations 

H = H'H = -{1-wb)pb-{^?, U'U = H'H. (42) 

The differentiation gives the autonomous system of 
equations 



x' = — 3a; + ^ 2;z + I X c(a:, y). 



\xz + ^yc{x,y), 



(43) 



{l + v^)z' = SX^xy ~ ^ zc{x,y), 
v' = ^vc{x,y), 



where c{x, y) = {1 + 'Wb){1 — x^ — y'^) + 2x^ . The quantity 
z is strictly constrained by the condition 



^ /I I 2n 2 Vo 2 



(44) 



which is the direct consequence of hyper-trigonometry: 
cosh^ q — sinh^ 9=1- This constraint makes system P5)) 
overdefined, since z is completely given by y, w and pa- 
rameter = Vo/Vq. Nevertheless, the system allows us 
to get a complete analysis of critical points in the simplest 
manner. What is of our interest? It is the projection of 
trajectory in the {x,y,v} 3D-space to the 2D-plane of 
{x,y}- 



A. Late times 

At present, the cosmological constant makes a sig- 
nificant contribution to the Hubble constant, i.e. v — 
Uq/U ~ 1. At late times of evolution just before the 
present, we put v ^ 1. Then, z can be excluded by 



(45) 



This limit means that the cosmological constant can be 
neglected, while the field has a large value, so that the 
hyper-sine can be approximated by a single exponent. 
Therefore, at late times we arrive to the analysis of ex- 
ponential potentials given in Indeed, under ([^5)1 
system (^5)) is reduced to the system for the exponential 
potential. The analysis of critical point in [llj gave the 
following physically meaningful properties: irrelevant of 
normalization of potential Vq there are stable scaling at- 
tractors in the plane of {x, y}; these attractors appear at 
> 3(1 + wb), so that = AV3(1 + wb) > I 
and W0 ~ Wb- The attractor is the stable node at 
A^ < 24(1 + wbY /{I + 9ws), otherwise it is the stable 
spiral focus. The position of attractors are given by 



3 1 -f Wb 



3(1-^1) 1 



(46) 



(41) that fixes z according to ([45)1 . i.e. Zc = 



Thus, at late times just before the present, the 
quintessence follows the scaling behavior independently 
of its initial conditions. 



B. Future 

Since function c{x,y) takes positive values at a; 7^ 0, 
?/ ^ 0, i.e. at presence of scalar field, quantity v grows 
in accordance with its differential equation in ()43|) . of 
course. Hence, in future we get v ^ 1. Then, yields 



z ^ Zi, 



Vo 



(47) 



i.e., z is frozen at z = z^,. Therefore, we get the system 
for the plane {x,y} in (|43p with z^, being the external 
parameter. The critical points are posed at the following 
sets: 
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I. Scalar field is absent, 

(i): Xi, = 0, 



2/* = 0, 



(48) 



so that the linearized equations in vicinity of critical 



point, i.e. with x 
system 



with the matrix 



B 



y* 



y, result in the 
(49) 



B 



having eigenvalues 



3{wb - 1) 2* 
-z* 3{wb + 1) 



(50) 



vi = 



V2 



i(3wB 



2) 



(51) 



that implies the critical point corresponds to instability 
due to ^tV2 > at wb > 0. Anyway, critical point (l48l) 
is a saddle at < 3 and wb = 0, which is of practical 
interest to look at the present time and future universe. 
We notice that actually the baryonic matter has a small 
pressure, which can be neglected in the universe evolu- 
tion, i.e. Wb +0. At \z^,\ > 3 eigenvalues (fSTI) satisfy 



J 



1^1 = 1^2' ^'^'^ have unstable focus at wb > or center 
at Wb = 0, while at wb < the focus becomes stable. 

II. Critical points of general position are given by 



(ii): 



(52) 



with additional symmetry over the following permuta- 
tions: A !->■ {Xi, ^ -Xi, and ~y*}, B ^ {x^, ^ y*}, 
C 1-^ {the product of operations A and B}. So, taking into 
account the symmetry, ([5^ covers 4 related sets. The ac- 
tion of permutation A conserves the stability properties, 
while action of B changes them. This fact becomes clear, 
for instance, at < 3. Indeed, altering the sign for 
X implies the interchange of two equivalent branches in 
the potential according to <-^- —(/) or formal altering the 
sign of quantity U, both of which change the sign of y, 
too. So, action A cannot influence the stability prop- 
erties. The action of B interchanges fractions of kinetic 
and potential energies in the energy density, that can be 
physically essential. 

The analysis of linear stability involves matrices for 
([5^ including action of symmetry A and for ([5^ con- 
verted by action of symmetry B, with corresponding 
eigenvalues {v}, 



^ I f WBiV^~zl-3)-2,/9~zl 
^ V -Zi,WB 



-Z^WB 



-wb{V9^ + 3) - 2^9^ 



■4 = ~iWB - V9^> 



^2 - -v^ 



(53) 



Bb = 



-u;b(V9+^+3) + 2V9^ 



-Zi,WB 



-Zi,WB 



Wb 



(V9^-3) + 2y9^ 



-3wb -I- ^9 - zl, 



(54) 



For real critical points, i.e. at |z^| < 3, and at wb ^ 0, set 
([5^ is stable node (both eigenvalues are negative) , while 
set (|52l) affected by action B gives saddle or unstable node 
depending on the sign of eigenvalue i.e. a balance 
between values of iwB and ^JW 



At \Zi,\ > 3 critical points ([5^ take complex values, 
while matrices satisfy Ba = B^, i.e. they are complex 
conjugate to each other as well as Xi, — y*. The same is 
valid for eigenvectors: 



3- 



z* 



3+ 



while e\ 2 — {^12}*- Basis (|55p is orthonormal: 



(55) 



„a a _ r 



The linearized system has solutions 



Xit 

y* 



(56) 



where Uj stand for some initial data. Solutions (|56p are 
complex- valued. By our investigation, (j56p are irrelevant 
to physical quantities in question at \zi,\ > 3. 

III. Tuned scaling appears at the special value of pa- 
rameter 2 



Wg. Then, there is the critical line 



(iii): = X* 



1 - Wb 



1 



WB 



(57) 



which is X-invariant, i.e. independent of slope in the 
potential. It is spectacular that tuned point given by 
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Zt = Zc, x^, = Xc, Hi, = He corresponds to the exact 
scaling solution found in Section II. 

At (|57p the linear analysis of perturbations gives two 
eigenvalues: vi — and 1/2 — '2wb, so that zero eigen- 
value corresponds to the line itself, while the positive one 
indicates instability at w_b > in vicinity of the line as a 
whole. This fact is in agreement with the study of criti- 
cal points in the case II, since (|57p contains the point of 
(|52| at < as well as it does the point at u;_b > af- 
ter the action of symmetry B, so we correspondingly get 
eigenvalues (j53p and (|54p at z^, z^ yielding the same 
result above. Therefore, at wb > the scaling solution of 
section II is unstable in future, while it does at wb < 0. 

IV. The boundary circle 



(iv): x'^ + 



1 



(58) 



is conserved by the autonomous system. This fact implies 
that at |z^| > 3 and wb > 0, when there are no stable 
critical points, the system approaches the boundary circle 
in future. 

Summarizing the results, we have found that in the 
presence of cosmological constant the scalar field with 
the hyper-cosine potential approaches the stable attrac- 
tor at late times just before the present day (the time of 
cosmological constant becomes visible) at the potential 
slope > 3(1-1- Wb), so that the energy of field scales 
like that of matter with the state parameter wb, while 
the fraction of field energy depends on the slope value. 
However, this attractor exhibits the strange behavior in 
future, i.e. under the dominance of cosmological term. 
Then, the balance of remnant energy depends on the pa- 
rameter Zi, determined by two quantities: i) the weight of 
potential normalization with respect to the energy den- 
sity due to the cosmological constant and ii) the slope, 
in accordance with (|T7)l . So, at < 3 and wb ^ the 
dynamical part of scalar field energy reaches the other 
scaling attractor, which dominates over the matter in 
accordance with ([55]) (see Fig. [T^) and gets the state 
parameter 



1 - 
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(59) 



At \zi,\ < 3 and wb < we get an interplay of two 
attractors with ((5T|) and ((53| : at 3wb + ^79 — < 
the scalar field relaxes at the minimum of kinetic and 
potential energy (|48p (see Fig. [T|d), while, otherwise, at 



3wb + — z1 > the scalar field dominates over the 
matter at the same point of ([53|) and (l59|) (see Fig. [TJ:). 

At > 3 we get three cases: i) at wb > the 
quintessence in future vibrates at the boundary circle as 
the limit cycle (see Fig. Wi;), ii) at wb = the field cycled 
around the center being the point of minimal kinetic and 
potential energy (see Fig. [If), ui) a.t wb < the field 
relaxes to the minimum (see Fig. HJi). 

The evolution kinds of scalar field in the plane of {x, y} 
are illustrated in Fig. [T]at A = 20 and various sets of z^, 
and Wb- at z* = 2.9 we put wb = 0.2 in Fig. [T]a), wb = 



—0.5 in b), and wb = —0.2 in c), while at = 10 we 
put Wb — —0.2 in Fig. [T]d), wb = 0.2 in e), and wb = 
in f). At all of tries the evolution starts at xq — —0.8, 
yo = 0.4, and trajectories move clockwise. We certainly 
see that trajectories approaches the late time attractor 
at appropriate {xc,yc}, which are negative in all cases 
except b), when they are positive. The future attractors 
in Fig. [U a) and c) are posed at {a;*,?/*} with opposite, 
negative, sign of (|52|) . The attractors at b) and d) stand 
in the minimum, while limit cycles of e) and f) are at the 
border and around the minimum, correspondingly. 

Thus, in the linear analysis we have classified the late 
time and future attractors for the quintessence with the 
specified kind of potential relevant to the case of nonzero 
cosmological constant. 

However, this analysis falls in special degenerate case 
of z* = 3, wb = 0. 



C. Degenerate case 

At Z-^ — wb — and v ^ 1 the analysis of future 
evolution becomes nonlinear, since, after the transforma- 
tion to variables a — x'^ + y'^ and t = x/y, autonomous 
system (^5)) is reduced to 



, , , 1 - 

cr = 3(7 cr- 1) — 2, 

1 -|- 



(60) 



which can be solved explicitly. Indeed, the integration 
for T results in 



(61) 



where 01 = 3(A^— A'o), and Aq corresponds to some initial 
data. Hence, r ^ 1 at the end of evolution, i.e. at 
N +00. Then, 



nnl^^ = hr^^+^,ln(^-01,) 



(62) 



where cro < Ij = ^(1 + i)- Therefore, at A^ ^ -t-00 
we find 

1 — (T (Tn 

In ^ -21nAf -00, 

1 - CTO cr 

that implies 

cr ^ 1, (63) 
and the attractor is posed at the boundary circle, 

(64) 



1 



Xi, = y^, = —= or x^, ^ y^, 
v2 



1 



If the initial data CTq = 1, the quantity a does not 
evolve, while r approaches the attractor r^, = 1. 
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In addition, we have explicitly found that the second 
critical point cr^ ~ 0, i.e. Xi, = y-^ = is unstable. 

The character of attractor ([64|) is illustrated in Fig. [21 
where trajectories move clockwise. 



1 








-0.75 - 0.5 - 0.25 0.25 0.5 0.75 1 



FIG. 2: The attraction of trajectories to the critical points at 
the boundary circle of {x, j/} plane in the degenerate case. 

Thus, attractors (|64p exhibit the stable behavior in ap- 
propriate semicircles, while the line connecting the criti- 
cal points has the instability versus perturbations. 



arbitrary to enlarge the mass up to reasonable values in 
the physics of Standard model. However, huge values of 
Zi, involve extremely frequent oscillations of quintessence 
in the nearest future, that is in contradiction with the 
present smooth evolution of universe. Therefore, wc ex- 
pect that a viable model includes the quintessence mass 
of the order of Hubble constant today. 

B. Restriction to the slope 

The late time scaling of quintessence results in a fixed 
fraction of quintessence energy in the budget of universe 
with respect to other matter irrespective of the evolu- 
tion stage: the dust or radiation fix close values of frac- 
tions. However, the fraction of nonbaryonic matter is 
constrained due to measured and primordial abundances 
of light elements caused by Big Bang Nucleosynthesis 
[Tol . Then, the slope of potential should be quite 

large to suppress fi^ ^ 0.13, so according [ll| one gets 

> 20. (68) 

Next, a role of quintessence field during inflation 
actually was analyzed in pT| , since the hyper-sine in fact 
coincides with the exponential potential at large values 
of field. The problem is a relic abundance of quintessence 
after inflation, that should be small in order to conserve 
the standard scenario of nucleosynthesis. Appropriate 
restrictions in various schemes of inflation are given in 




IV. PHENOMENOLOGICAL POINTS 
A. The mass 

The potential of quintessence suggests the mass 

m^ = ^ =47rGA2yo. (65) 
d(j)^ 0=0 

For the exact scaling solution in Section II we get 

ml^^-{l^wl)Hln^, (66) 

and at = and ~ 0.7 of practice the mass is 
determined by the current value of Hubble constant, i.e. 
it is extremely small as well as the energy scale of cos- 
mological constant, that is beyond a natural reason. Al- 
though, such the mass could argue for both the present 
acceleration in the universe expansion and the scale of 
cosmological constant. 
Generically, we get 

2tt 1 
= — VoG zl ^ - H^Qa zl (67) 

So, the mass of scalar quintessence scales as the present 
day Hubble constant with the factor of z^,, which could be 



C. Initial conditions 

Attractors mean a slow dependence of late time evolu- 
tion on initial data for the quintessence. The character 
of regulation is illustrated in Fig. [3l The set of tries 
exhibits the following general features: 

• At small initial fraction of quintessence energy, it 
is frozen to a moment, when it approaches an ap- 
propriate scaling value in order to start the tracker 
behavior at late times. 

• At large initial fraction of quintessence energy, it 
rapidly falls in order to frozen and wait for a mo- 
ment of tracker way at late times. 

• In future, vibrations of quintessence at z^, > 3 and 
wb ^ determine an average value of dynamical 
parameter for the equation of state {w^j,), which is 
independent of initial data, whereas — 1 < (w^) < 
wb at wb > or {w^) = at wb = (see Fig. 
[3] a, b, c, d). At Wb < vibrations determine the 
effective {w^) > wb (see Fig. [Hlg, h). 

Further observations repeat general properties of fu- 
ture attractors. 
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FIG. 3: Approaching the scaling under different initial data: the logarithm of energy In p in arbitrary units versus the e- folding 
TV = Ina(i) with a casual shift of present day scale factor a{to); the quintessence and baryotropic matter are represented by 
thick solid and dashed lines, correspondingly, while the dynamical part of quintessence without the cosmological term is given 
by thin solid line. Initial data for the suppressed and dominated quintessence are respectively shown in left and right panels: 
a) and b) at A = 20, wb = and 2;* — 50; c) and d) at A = 20, wb = 0.2 and = 10; e) and f) at A = 20, wb ~ 0.2 and 
= 2.9; g) and h) at A = 20, wb = -0.2 and z, = 10. 
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FIG. 4; The state parameter of quintessence versus the 
e-folding A'' of evolution scale: changing the stable value of 
late times at = at = 50 and A = 20. 



D. Equation of state 

At late times, the attractor causes the ratio of 
quintessence pressure to the energy density Wcj, is sta- 
bilized infinitely close to the value of parameter wb for 
the matter. However, the situation changes, when the 
cosmological term comes to dominate, and moves to 
— 1. An example of relaxation is shown in Fig. |4]for 
the quintessence vibrating around the minimum point 
(see Figs. [T]f and [3] a, b). The present day e-folding is 
arbitrary in Fig. [51 The magnitude of deviation from 
the limit of —1 and vibration period depend on the po- 
tential parameters. The picture analogous to Fig. [Hwas 
observed in \2§\ with similar values of parameters. 

It is clear that vibrations are absent if < 3. Any- 
way, the parameter of quintessence state u;^ rapidly ap- 
proaches the vacuum value. Nevertheless, it would be 
interesting to see a relative significance of quintessence 
with respect to the matter. So, general consideration of 
attractors and Fig. [3] demonstrate that the quintessence 
fraction can dominate or be suppressed depending on the 
potential parameters, the value of z^,. It is evident that 
in the case of reaching the boundary circle in the {x, y} 
plane the effective, average value for the state parameter 



is (w^) = 0, while the same value is clearly observed also 
in the case of wb = and \Zi,\ > 3 (see Fig. |3]a and b). 

Thus, we get the definite understanding of phenomeno- 
logical properties for the evolution of quintessence with 
the specified kind of potential in the presence of cosmo- 
logical term. 



V. CONCLUSION 

In this paper we have found the potential of scalar 
field quintessence, which gives the exact solution for the 
scaling evolution of flat universe in the presence of cos- 
mological constant. The scaling behavior is consistent 
with the current empirical observations. 

We have investigated the stability of scaling behavior 
versus the variations in the slope and normalization of 
potential as well as in initial data. The analysis has re- 
vealed two kinds of attractors. The late time attractor 
just before the cosmological constant is coming to play, 
is independent of normalization, and it is determined by 
the slope, that is consistent with the well-known result 
for the exponential potentials [11], representing the limit 
of large field for the potential found in the paper. The fu- 
ture behavior of quintessence under the dominance of cos- 
mological constant depends on both the ratio of potential 
normalization to the vacuum energy density and slope in 
the special combination denoted by parameter Zi,. Gener- 
ically, the future attractor differs from that of late time. 
So, the late time attractor reveals the strange behavior. 
We have classified the future attractors by their charac- 
ter and stability in linear analysis. The degenerate case 
of nonlinear dependence has been solved explicitly. Some 
phenomenological items have been considered, too. 

We conclude that analysis of scaling attractors can be 
useful for classifying the quintessence behavior at late 
times and in future. 

This work is partially supported by the Russian Foun- 
dation for Basic Research, grant 04-02-17530. 
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